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ABSTRACT 

We present Chandra studies of the X-ray binary (XRB) populations in the bulge and ring regions of the 
ring galaxy NGC 1291. We detect 169 X-ray point sources in the galaxy, 75 in the bulge and 71 in the 
ring, utilizing the four available Chandra observations totaling an effective exposure of 179 ks. We report 
photometric properties of these sources in a point-source catalog. There are sa 40% of the bulge sources and 
w 25% of the ring sources showing > 3(7 long-term variability in their X-ray count rate. The X-ray colors 
suggest that a significant fraction of the bulge (w 75%) and ring (ss 65%) sources are likely low-mass X-ray 
binaries (LMXBs). The spectra of the nuclear source indicate that it is a low-luminosity AGN with moderate 
obscuration; spectral variability is observed between individual observations. We construct 0.3-8.0 keV X- 
ray luminosity functions (XLFs) for the bulge and ring XRB populations, taking into account the detection 
incompleteness and background AGN contamination. We reach 90% completeness limits of « 1.5 x 10^^ and 
w 2.2 X 10^^ erg s"' for the bulge and ring populations, respectively. Both XLFs can be fit with a broken 
power-law model, and the shapes are consistent with those expected for populations dominated by LMXBs. 
We perform detailed population synthesis modeling of the XRB populations in NGC 1291, which suggests 
that the observed combined XLF is dominated by an old LMXB population. We compare the bulge and ring 
XRB populations, and argue that the ring XRBs are associated with a younger stellar population than the bulge 
sources, based on the relative overdensity of X-ray sources in the ring, the generally harder X-ray color of the 
ring sources, the overabundance of luminous sources in the combined XLF, and the flatter shape of the ring 
XLF 

Subject headings: galaxies: individual (NGC 1291) — galaxies: luminosity function — X-rays: binaries — 
X-rays: galaxies 
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1. INTRODUCTION 

Galactic rings are a peculiar phenomenon among galaxies. 
A prototype ring galaxy is the famous Cartwheel galaxy, with 
a spectacular ring luminous in the infrared (IR), optical, UV, 
and X-ray bands. The formation of ring-like structures in 
galaxies is not fully understood. The two most common theo- 
ries are the collisional origin and the resonance origin. A col- 
lisional ring could be created by a head-on collision between 
a compact companion galaxy and a larger disk system; a col- 
lisional ring galaxy usually contains an off-centered nucleus 
and an asymmetric ring structure, and there is a significant 
chanc e of finding a companion galaxy close to the ring galaxy 
(see Apple ton & Struck-MarceUll 19961 for a review). A reso- 
nance ring could be produced by the outer Lindblad resonance 
of spiral bars; such a system is preferentially observed in 
barr ed galaxies, and the res ulting ring is generally symmetric 
(see lButa & Combelll 19961 for a review). There are a few sys- 
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tems that are difficult to explain w ith these two theories, such 
as the Hoag's Object (Hoag 1950), for whic h alternative for- 
matio n scenarios have been proposed (e.g., iFinkelman et alj 
l20Tlh . 

Galactic rings are expected to exhibit significant star for- 
mation compared to their progenitor galaxies mainly due to 
the enhanced gas densities in the rings, regardless of the 
formation mechanism. Surveys of ring galaxies in optical. 
Ha, IR, and radio have revealed that th e rings are generally 
rich in Hi g as and Hll regions (e.g., van Driel et al. '19881; 
Marston & A ppleton 1995;"Phillips 1996; Ryder etal. 1991 
Grouchy et alJ I2010). With Chandra's subarcsecond angular 



resolution, it is now possible to study in detail the X-ray bi- 
nary (XRB) populations in nearby galaxies, e.g., low-mass 
X-ray binaries (LMXBs) that are associated with old stellar 
populations and high-mass X-ray binaries (HMXB s) that are 
associated with young stellar populations (see e.g., Fabbiano| 
I2OO6 for a review). Therefore, X-ray observations have be- 
come an independent approach for constraining the nature of 
the underlying stellar populations by means of X-ray lumi- 
nosity functions (XLFs) and X-ray colors of the XRBs. How- 
ever, little is known regarding the X-ray properties of galac- 
tic rings, except for the Cartwheel galaxy, which has an ex- 



X-ray sources (ULXs; Lx > 10-^ erg 
2004) in the star-formin 
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NGC 12 9 1 is a remarkable ring galaxy, first noted 
by iPerrind (Il922h. It was c lassified as (R)SB(s)0/a 

i RC3; Ide Vaucouleurs et al.l Il991h . and is virtually face-on 
MeboldetalJll979h . The images of NGC 1291 in the UV, 
optical, and near-IR (NIR) are shown in Figure [l] The 
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galaxy includes a fainter outer ring and a prominent bulge, 
which was considered as a p rototypical pseudobulge (e.g., 
iKormendv & Kenniculj 120041). At a d istance of 8.9 Mpc, 
NGC 1291 has Mb,o=-20.6 ddeVaucouleurs 197 5]). The D25 
ellipse of NGC 1291 has a size of 9.8' x 8.1' (iKennicutt et al.l 
12003) . and the nuclear position of the galaxy is ajaooo = 
03''17™18.6% 5j2ooo.o = -41°06'29" (lEvans et al.l l2010). The 
Galactic column density alon g the line o f sight to NGC 1291 
is A^H = 2.1 X 10^" cm-2 ( Dick ey & Lock man 1990)."' 

The ring of NGC 1291 does not appear as energetic as that 
of the Cartwheel galaxy. There were only two ULXs re- 
porte d in the ring based on the first Chandra observation in 
2000 ( Swartz et al]l2004 . The total st ar formation rate (SFR ) 
of NGC 1291 is only 0.4 Mq yr"' jKennicutt et al.ll2003h . 
compared to that of 67.5 Mq yr~' for the Cartwheel galaxy 
(|Higdon 1993). Given its barred spiral nature, the relatively 
symmetric structure of the ring, and the absence of any ap- 
parent companion, NGC 1291 is likely a resonance ring sys- 
tem, while the C artwheel galaxy is most likely a c ollisional 
ring system (e.g.. Appleton & Stru ck-Marcelllll996l) . Intense 
star formation still appears present in the ring of NGC 1291, 
which is obvious from the UV image (Figure[T]l. On the con- 
trary, the bulge is much redd er than the ring, and d ominates 
in the NIR image (Figure [l). Ivan Driel et al.l d 19881) mapped 
the Hi distribution in NGC 1291, and revealed that Hi gas is 
concentrated in the ring region. The ring also contains dis- 
crete Hl^j;egions, while the bulge has filaments of ionized 
gas (ICrocker et al.lll996l) . Given the UV/optical/radio data, 
the rin ^ appears t o host a younger stellar population than the 
bulge. 'Nol l et akl (l2009h studied the star-formation history of 
many galaxies including NGC 1291, using spectral energy 
distribution fitting techniques. The best-fit model for NGC 
1291 (Noll, private communication) includes two exponential 
starburst events: one happened 10 Gyr ago with an e-folding 
timescale of 0.25 Gyr and the other happened 200 Myr ago 
with an e-folding timescale of 50 Myr The stellar mass of the 
old population is « 330 times higher than that of the younger 
one. Although the fitting results were derived for the entire 
galaxy and may have large uncertainties, they are consistent 
with the UV/optical/radio data if the ring hosts the younger 
stellar population. 

NGC 1291 had been observed with Chandra twice in June 
and Nov ember 2000, with a total exposure of « 80 ks. A first 
study by llrwin et al.l (l2002h was limited to the bulge region. 
lKilgardetal.l(l2005h presented a point-source catalog based 
on these observations, and found that the XLF and X-ray col- 
ors suggest old XRB populations in the galaxy. This study 
however, did not discriminate between bulge and ring XRBs. 

We obtained two additional Chandra observations of NGC 
1291 in May 2010, with a total exposure of w 120 ks. The new 
observations cover the majority of the ring region, allowing 
us to investigate the XRB populations in the bulge and ring 
respectively. NGC 1291 is much closer than the Cartwheel 
galaxy (« 10 Mpc compared to « 120 Mpc), and thus we are 
able to probe the ring XRB population down to an unprece- 
dented luminosity limit (w lO-'^ erg s~'). We compare the ring 
XRB population to the bulge population and test whether the 
galactic ring hosts a younger stellar population from an X-ray 
point of view. The deep observations on the bulge also allow 
us to compare its LMXB population to those of the recently 
acquired deep exposures of elliptical galaxies, and establish if 
all LMXB populations are indeed compatible. 



This paper is organized as follows. We describe the Chan- 
dra observations of NGC 1291 and the data reduction proce- 
dure in Section 2. In Section 3 we construct an X-ray point- 
source catalog. X-ray source properties are presented in Sec- 
tion 4, including X-ray source radial profile, source variabil- 
ity, and X-ray colors. The nuclear source and ULX candidates 
are also discussed. In Section 5, we calculate and parameter- 
ize the XLFs for the bulge and ring regions, respectively, tak- 
ing into account the detection incompleteness and background 
AGN contamination. We model the observed XLF using pop- 
ulation synthesis (PS) models in Section 6. In Section 7, we 
discuss implications of the X-ray properties and we argue that 
the ring hosts a younger stellar population than the bulge. We 
summarize our results in Section 8. 

2. OBSERVATIONS AND DATA REDUCTION 

NGC 1291 has been observed with the S3 ch ip of Chandra 
Advanced CCD Imaging Spectrometer (ACIS: iGarmire et al.l 
2003) at four epochs. The four observations are listed in Ta- 
ble [T] along with their exposure times, ranging from 37 ks to 
70 ks. The focal-plane temperature was -120°C for all the ob- 
servations. We reduced and analyzed the observational data 
using mainly the Chandra Interactive Analysis of Observa- 
tions (CIAO) tools. " We used the CHANDRA_REPRO script 
to reprocess the data with the latest calibration. The back- 
ground light curve of each observation was then inspected and 
background flares were removed using the DEFLARE CIAO 
script. The flare-cleaned exposure times are also listed in Ta- 
ble [1] there was a significant flare in observation 2059 that 
reduces its usable exposure by « 40%. The total usable expo- 
sure for NGC 1291 is 178.6 ks. 

We registered the astrometry of the observations using the 
CIAO tool REPROJECT_ASPECT. We created a 0.3-8.0 keV 
image for each o bservation and se arched for sources using 
WAVDETECT (Freem an et al.ll2002l) at a false-positive proba- 
bility threshold of 10"^. These source lists were used to reg- 
ister the observations to the astrometric frame of observation 
11272 with the CIAO script REPROJECT_ASPECT, adopting 
a 3" matching radius and a residual rejection limit of 0.6". 
We reprojected the registered event lists to observation 1 1272 
using REPROJECT_EVENTS, and then merged all the obser- 
vations to create a master event file using DMMERGE. The 
ACIS-S3 chip has different pointings for the four observa- 
tions, and the average aim point (weighted by exposure time) 
is aj200().o = 03M7™09.3P, ,5j2ooo.() = -41°06'25.1". 

We constructed images from the merged event file using 
the standard ASCA grade set (ASCA grades 0, 2, 3, 4, 6) 
for six bands: 0.3-8.0 keV (full band), 0.3-2.0 keV (soft 
band), 0.5-2.0 keV (conventional soft band), 2.0-8.0 keV 
(hard band), 0.3-1.0 keV (soft band 1), and 1.0-2.0 keV (soft 
band 2). For each observation, we created exposure maps in 
these band s following the basic proce dure outlined in Sec- 
tion 3.2 of Hornschemeier et al] (12001b . which takes into ac- 
count the effects of vignetting, gaps between the CCDs, bad- 
column filtering, bad-pixel filtering, and the spatially depen- 
dent degradation in quantum efficiency due to contamina- 
tion on the ACIS optical-blocking filters. A photon index of 
r = 1 .7 was assumed in creating the exposure map, which 
is a typical value for photon indices of X -ray binaries (e.g., 
llrwin et ani2003l : iBrassington et al.ll2010 ). Merged exposure 
maps were then created from the exposure maps of the indi- 
vidual observations. 



See http://heasarc.nasa.gov/cgi-bin/Tools/w3nh/w3nh.pl. 



See |http://cxc.harvard.edu/ciao7] for details on CIAO. 



Fig. 1 .— Images of NGC 1291 in the GALEX near-UV (NUV, Aeff = 2267 
image is 14' X 14', centered on the galactic nucleus. The ellipses outline the 
dominated by young stars. 
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Fig. 2. — Chandra "false-color" image of NGC 1291. This image is acolor 
composite of the exposure-corrected adaptively smoothed images in the 0.3- 
1.0 keV (red), 1.0-2.0 keV (green), and 2.0-8.0 keV (blue) bands. The cross 
symbol shows the center of the galaxy. The white dashed ellipse shows the 
D25 region I Kennicutt et al. 2003). The cyan and magenta ellipses outline the 
bulge and ring regions of the galaxy, respectively. 



We created adaptively smoothed images from the raw im- 
ages using the CIAO tool CSMOOTH. Exposure-corrected 
s moothed images were th en constructed following Section 3.3 
of iBaganoff et al.l (l2003h . We show in Figure |2] a color com- 
posite of the exposure-corrected smoothed images in the 0.3- 
1.0, 1.0-2.0, and 2.0-8.0 keV bands. 

3. POINT-SOURCE CATALOG 



A), Digitized Sky Survey red (RG610), and 2MASS K bands. The size of each 
bulge and ring regions (see Section 3 for details). The UV ring appears to be 

We ran WAVDETECT on the full-band image to search for 
X-ray point sources, with a "\/2 sequence" of wavelet scales 
(i.e., 1, ^/2, 2, 2^2, 4, 4^2, 8, 8%/2, and 16 pixels) and a 
false-positive probability threshold of lO"*". At this threshold, 
we expect approx imately two false detections given the size 
of the image (e.g.. lKim et"ani2004l) . 

There are 169 X-ray sources det ected in the full band. 
We utilized the ACIS EXTRACT (AE: iBroos et al.l l2010') pro- 
gram to derive the source photometry. The main advantage 
of AE compared to a simple circular-aperture extraction ap- 
proach is that AE uses a polygonal extraction region to rep- 
resent the point spread function (PSF) for each observation 
of each source individually, which is more suitable for mul- 
tiple observations with different aim points and roll angles. 
Briefly, AE constructs source extraction regions that approxi- 
mate w 90% encircled-energy fraction (EEF) contours of the 
PSFs at 1 .50 keV for each exposure. Smaller extraction re- 
gions (as low as « 60%) were used for sources in crowded 
regions. Background counts were extracted in an annular re- 
gion around each source excluding any pixels that belong to 
detected sources; the number of extracted background counts 
was required to be at least 100 to reduce the uncertainty from 
the background estimation. Encircled-energy fractions at sev- 
eral energies (ranging from 0.28 to 8.60 keV) were provided 
by AE, and thus we derived aperture corrections for all the six 
X-ray bands listed in Section |2] via interpolation. Aperture- 
corrected net (backgroun d-subtracted) source counts and the 
corresponding \a errors ( Gehrelsl ll986') for these bands were 
then computed; the errors were propagated through the errors 
of the extracted sources and background c ounts fo llowing the 
numerical method described in §1.7.3 of iLvonsI ( 11991) . As 
the source detection was performed only in the full band, we 
consider a source as undetected in a given band (other than the 
full band) if it is less significant than 3ct, when calculated us- 
ing the net counts and Ict error in that band. An upper limit on 
the source counts was der ived for an undete cted source, using 
the Bayesian approach of iKraft et"an (1199 Ih for the 99% con- 

Source searching was performed only in the full band as we focused on 
the study of the full-band XLFs. 
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fidence level if the extracted number of counts is less than 10 
in the given band; otherwise, \ ye calculated a 3a upper limit 
following the Poisson statistic (lGehrelslll986l) . 

An X-ray spectrum for each source was extracted by AE 
along with the spectrum of its local background. We per- 
formed sgectral fitting of each source using XSPEC (Version 
12.6 0:lArnaudlll996t) with the Cash fitting statistic (CSTAT; 
ICashlll979h . An absorbed power-law model (TBABS*POW) 
was employed for sources with more than 30 full-band counts; 
for the other sources, the spectra were modeled with a simple 
power law with the photon index fixed at F = 1.7. The ob- 
served full-band fluxes were derived from the best-fit spectra, 
which were then converted to full-band luminosities. We ap- 
plied a smaller correction factor (1.05) to the observed full- 
band luminosities to account for the Galactic absorption. 

Of the 169 point sources, 75 are in the bulge region and 
71 in the ring area. We defined the bulge and ring regions 
based on the UV and NIR images shown in Figure 1. We 
examined the surface brightness contours of the /iT-band im- 
age and selected a contour that encloses w 80% of the light 
within the D25 ellipse. This contour can be well approxi- 
mated by an 4.7' x 3.7' ellipse centered on the galactic nu- 
cleus, which was used to defined the bulge boundary. The ring 
region is defined by two ellipses centered on the galactic nu- 
cleus. We fixed the position angle (70°) and axis ratio (1.15) 
o f the ellipses according to the morphology information given 
in lde Vaucouleursl(ll975l) and only allowed the sizes of the el- 
lipses to vary so that the ring region includes the majority of 
the prominent UV ring features. We note that 85% of ring 
region is covered by the Chandra exposure (see Figure|2]i. 

The spatial distribution of the X-ray sources is displayed in 
Figure [3^, and their luminosity distribution is show in Fig- 
ure |3j'. The full-band luminosities range from 2.7 x 10^^'' to 
2.3 X 10^' erg s~'. We listed basic photometric properties for 
the 169 sources in Table 2, with details of the columns given 
below. 

1 . Column 1 : the source number, listed in order of increas- 
ing right ascension. 

2. Columns 2 and 3: the right ascension and declination 
of the X-ray source, respectively. 

3. Column 4: the off-axis angle of the X-ray source in ar- 
cminutes. This is the distance between the source posi- 
tion given in columns 2 and 3 and the average aim point 
(see Section|2|l. 

4. Column 5: the radial distance of the X-ray source to the 
nucleus of the galaxy in arcminutes. 

5. Columns 6-23: the source counts and the corre- 
sponding Icr statistical errors or the upper limits 
on source counts for the six bands, 0.3-8.0 keV, 
0.3-2.0 keV, 0.5-2.0 keV, 2.0-8.0 keV, 0.3-1.0 keV, 
and 1.0-2.0 keV, respectively. Upper limits are indi- 
cated as a "-1 .0" in the error columns. 

6. Column 24: the full-band effective exposure time de- 
termined from the exposure map, which has been cor- 
rected for vignetting and quantum-efficiency degrada- 
tion. 

7. Columns 25-30: the soft and hard X-ray colors and the 
corresponding Icr statistical errors; see Section |43] be- 
low for details. 



8. Column 31: the photon index derived from spectral fit- 
ting, or set to F = 1 .7 if the number of full-band counts 
is smaller than 30. 

9. Column 32: observed full-band X-ray luminosity, in 
units of erg s~' . 

10. Column 33: the statistical significance of the full-band 
flux variability; see Section l43] below for details. It is 
set to "-1 .0" if the flux variability is not available. 

11. Column 34: region flag for the source. The flag is 
"Bulge" if the source is located in the bulge, "Ring" 
if the source is in the ring, or "Other" if the source is 
neither in the bulge nor in the ring. 

4. POINT-SOURCE PROPERTIES 

4.1. Radial Profile 

We constructed a radial profile of the X-ray source surface 
density. Given the non-circular shapes of the ring and bulge, 
we calculated X-ray source densities in four elliptical annular 
regions in the bulge, four elliptical annular regions in the ring, 
and one region in between the bulge and the ring. The ellip- 
tical annular regions in the ring/bulge have the same position 
angle and axis ratio as the ring/bulge and are evenly spaced 
in their distances to the galactic nucleus. We corrected the 
number of X-ray sources in each region for the detection in- 
completeness and background AGN contamination (see Sec- 
tion |5T| below; the incompleteness correction was derived in 
the corresponding region). For this purpose, we only included 
sources that are outside the central 10"-radius region and have 
full-band luminosity Lx > 1 x 10^^ erg s"^ The corrected 
number of sources was then divided by the region area to ob- 
tain the surface density. The X-ray source radial profile is 
presented in Figure |4] There appears to be an overdensity 
of X-ray sources in the ring of the galaxy compared to the 
region in between the bulge and the ring. We also included 
in this plot surface brightness radial profiles for the GALEX 
NUV and 2MASS ^T-band data, calculated in the same re- 
gions as above. The surface brightness profiles are arbitrarily 
normalized so that the highest brightness has the same value 
as the highest source density. The NUV surface brightness 
is expected to trace the young stellar populations, while the 
^T-band emission is sensitive to the old stellar populations. In 
the ring region, the NUV and /T-band surface brightness pro- 
files differ significantly. The X-ray source radial profile in the 
ring does not continue with the trend shown in the bulge and 
the region in between, and it mainly follows the NUV surface 
brightness distribution, suggesting that the X-ray sources in 
the bulge and ring regions are associated with different stellar 
populations and the ring stellar population is likely younger. 

4.2. Central Point Source and Ultraluminous X-Ray Sources 

The nuclear point source (XID 84) of NGC 1291 has a po- 
sition of aj2ooo.o = 03^'17™18.58^ (5j2()oo.o = -41°06'28.7", and 
an observed full-band luminosity of 2.3 x 10^^ erg s"'. It has 
been classified as a low-ionization n uclea r emission-line re- 
gion (LINER; e.g.. lSmith etallllOOTh . and llrwin et all (|2002|) 
suggested that it is an obscured low-luminosity active galac- 
tic nucleus (AGN) based on the X-ray spectral analysis of the 
first two Chandra observations. We show the X-ray spectra 
of the source for the four individual observations in Figure |5] 
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Fig. 3. — (a) Point-source spatial distributions for tlie 75 bulge sources (filled dots), 71 ring sources (open dots), and 23 other sources (triangles). The dashed 
boundary lines shows the ACIS-S3 coverage of NGC 1291. (b) Full-band luminosity vs. distance to the nucleus of the galaxy. 
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Fig. 4. — Radial profile of the X-ray source surface density. The densi- 
ties were computed in four elliptical annular regions in the bulge (green filled 
circles), four elliptical annular regions in the ring (magenta open circles), 
and one region in between the bulge and the ring (black triangle). The x-axis 
shows the average distance of each region to the galactic nucleus. The surface 
densities have been corrected for detection incompleteness and background 
AGN contamination. Only sources outside the central 10"-radius region and 
having Lx > lO'^ erg s"' are included in the density calculation. The errors 
of the surface densities were calculated based on the Poisson errors of the 
number of sources in each bin (Gehrels 1986). The red and blue curves dis- 
play the surface brightness radial profiles (scaled to the highest X-ray source 
density) for the GALEX NUV and 2MASS Tf-band data, respectively. There 
appears to be an overdensity of X-ray sources in the ring region. (A color 
version of this figure is available in the online journal.) 



which were extracted by AE with polygonal regions that ap- 
proximate the « 89% EEF PSFs in the 0.3-8.0 keV band. Lo- 
cal background spectra were also extracted by AE in the same 
annular regions as those used in the photometry extraction. 
There are clearly strong soft X-ray excesses in the spectra of 
all the observations, which are typical among low-luminosity 
AGNs and are considered to originate from hot gas in the 



galactic center (e.g..' Ptak et alJfT999 l). We thus fit the spectra 
using XSPEC with an absorbed power-law (AGN) plus ther- 
mal plasma (hot gas) model (WABSl*APECH-WABS2*POW). 
Note that there could be multiple temperatures present in the 
gas, however, the limited signal-to-noise ratio (S/N) of the 
data cannot place strong constraints on a multi-temperature 
plasma. The spectral slope and column density for the power- 
law component were free to vary during the fitting. For the 
thermal component, we set the abundance to be Solar, and 
then we determined its temperature (O.lSilijj'Qj keV) and ab- 



sorption column density (A^h.i = 7.6_o 7 x 10" cm~") by fitting 
the stacked spectrum of the four observations. In the spectral 
fitting of the individual observations, the temperature and col- 
umn density parameters were fixed, and only the normaliza- 
tion was allowed to vary.'^ The best-fit parameters are listed 
in Table |3j also shown are the observed luminosities of the 
AGN component. The la errors were derived by varying one 
interesting parameter The intrinsic 2.0-8.0 keV X-ray lumi- 
nosity is w (2-3) X lO-''^ erg s"' after absorption and aperture 
corrections, indicating its low-luminosity nature. The nuclear 
AGN has an intrinsic power-law photon index of w 2 and an 
absorption column density of w 2 x 10'^ whic h are typ- 
ical for an obscured AGN (e.g., [Turner et al.lll997l) . The 3cr 
confidence contours for the photon index and column density 
parameters are shown in Figure |6l there is a « 3(T spectral 
variation between observations 795 and 1 1272. The X-ray lu- 
minosities of these two observations also differ at the ss 3cr 
significance level (Table 1). It has been noted that long-term 
variability in luminosity and spectral parameter s is common 
among obscured AGNs (e.g.. lRisaliti et al.ll2002l) . 

Besides the nuclear source, there are two more sources 
(XIDs 55 and 57) with observed full-band luminosity in ex- 

The extraction regions are different between observations, correspond- 
ing to the fa 89% EEF PSFs of the nuclear point source in individual obser- 
vations. Therefore, the amount of diffused gas within the extraction region 
varies. 
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Fig . 5 . — X-ray spectra of the nuclear source overlaid with the best-fit mod- 
els for the four individual observations. The spectra were fit with an absorbed 
power law plus thermal plasma model (WABS 1 *APEC+WABS2*POW). The 
dashed curves show the individual model components. For the thermal com- 
ponent, the temperature (0.15 keV) and column density (7.6 X 10^' c m~^) 
were fixed at the values derived from the stacked spectrum. See Section l4!2l 
for more details and see Table|3]for the fit parameters. 

cess of 10^^ erg s~', in the regime of ULXs. Their in- 
trinsic full-band luminosities after absorption corrections are 
2.1 X 10^^ erg s~' and 2.0 x 10-''' erg s"', respectively. One 
other source (XID 26), after being corrected for intrinsic ab- 
sorption (Nh = 1 .4 X 10^' cm"^), has an intrinsic full-band lu- 
minosity of 1.1 X 10^^ erg s"^ which should also be consid- 
ered as a ULX candidate. XIDs 26 and 55 are ring sources, 
whereas XID 57 is a bulge source. The two ring ULX candi- 
dates have been reported in the ULX catalog of Swartz et al. 
([2004) based on the data of the first observation (observa- 
tion 795). The bulge ULX candidate is a transient source, 
not detected in the first two observations . These ULX can - 
didates were not detected in the study of llrwin et al.l (|2002|) . 
which utilized only the first observation and focused on the 
bulge region. There is also a non-negligible possibility that 
these sources are background AGNs. The expected numbers 
of background AGNs with ULX luminosities if placed at the 
distance of NGC 1291 are « 0.05 and w 0.3 for the bulge and 
ring regions, respectively (see Section 15.11 below for details 
of the AGN surface density). The numbers of ULXs in NGC 
1291 and the Cartwheel galaxy (« 3 and « 20) do not appear 
to be linearly correlated with their SFRs (0.4 M© yr"' and 
67.5 Mq yr"'). Previous ULX surveys have found that the 
numbers of ULXs at low SFRs (< 0.4 Mq yr"') are signif- 
icantly higher than those extrapolated linearly from the high 
SFR end, suggestive of an a dditional popul ation of ULXs as- 
sociated with LMXBs (e.g.. lLiu et al.ll2006h . 

4.3. Source Variability 

X-ray flux/spectral variability is a common feature among 
point sources in galaxies, which is generally attributed 
to the ch ange of p h ysical p roperties of the accretion 
disks (e.g., Remillard & McClintocik '2006; Done et al.ir2007l; 
iBrassington et al. 2010.; Fabbiano et al. 2010Ik The four 
Chandra observations of NGC 1291 span « 10 years, al- 
lowing us to study the long-term flux variability of the point 
sources. 

For each source in each observation, we calculated its 
aperture-corrected full-band count rate (CR) based on the pho- 
tometry provided by AE; an upper limit on the count rate was 



Obs 1 1 273 




1.0 1.5 2.0 2.5 3.0 3.5 

Fig. 6. — 3(T confidence contours for the photon index and absorption 
column density parameters of the nuclear AGN in the three individual ob- 
servations, observations 795 (blue), 11272 (green), and 11273 (red). The fit 
parameters in observation 2059 have large uncertainties and is not shown in 
the plot for clarity. The crosses show the best-fit A'H.a and F values in obser- 
vations 795, 11272, and 11273, respectively. The larger confidence contour 
of observation 11273 than that of observation 795 despite of its longer expo- 
sure is due t o the lower flux/luminosity of the nucleus during this observation. 
See Section 1431 for more details. (A color version of this figure is available 
in the online journal.) 

calculated if the source was not detected in the observation 
(see Section[3]l. For a source that was detected in at least one 
individual observation and was covered by at least two obser- 
vations, we computed the maximum statistical significance of 
its fl ux variation between a ny two observations, defined as 
(e.g.. lBrassington et alJl2009HSeU et all201 ih : 

\CRi-CRi\ 

ayar = max,, , ( 1 ) 

/ 2 2 

where the subscripts / and j run over different observations, 
and (Tc«, acRj are the Icr errors of the count rates. We con- 
sider a source to be variable if its dyar parameter is greater than 
three (i.e., > 3(t variation). The flux variability for the point 
sources in NGC 1291 are shown in Figure [7] In the bulge 
and ring regions, 3^t})^% and 25^1^% of the sources are vari- 
able, respectively; the errors are the la Poisson uncertainties 
on the numbers of sources that are considered variable. We 
note that the central point source and the three ULX candi- 
dates are variable. There is no significant difference between 
the fraction of variable so urces in the bulge an d ring regions 
given the data available. iKilgard et al.l (l2005h reported that 
29% of the sources show long-term variability based on the 
first two Chandra observations, consistent with our findings 
here. The fraction of variable sources in NGC 1291 is in gen- 
eral agreement with those in some t ypical early-type galax- 
ies, e .g., NGC 3379 and NGC 4278 (IBrassington etalllEoOSl 
l2009h . 

For the luminous nuclear point source and the three ULX 
candidates, we also examined their spectral variability by fit- 
ting the spectra in individual observations. The spectral vari- 
ability of the nuclear source is discussed in Section l4~2l For 
the ULX candidates, their spectra were well fit by an absorbed 
power-law model. Although all these ULX candidates are 
variable in flux, none of them shows significant (> 3a) vari- 
ability in its spectral parameters. The two non-transient ULX 
candidates (XIDS 26 and 55) still have ULX luminosities in 
individual observations. 
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Fig. 7. — Significance of flux variability vs. observed full-band luminos- 
ity for the 50 bulge sources (filled dots), 28 ring sources (open dots), and 3 
other sources (triangles). Dotted line indicates the > 3<j significance level of 
the flux variability. Note that we calculated flux variability for only sources 
detected in at least one individual observation; there are insufficient data for 
variability study of the remaining sources. We studied spectral variability 
of the few ULX candidates, which show no significant variabiUty in spectral 
parameters (column density and power-law index). 

4.4. X-Ray Colors 

It is difficult to classify individual X-ray sources into differ- 
ent XRB types (e.g., LMXBs and HM XBs) given the limited 
X-ray data alone (e.g.. lFabbianoll2006l) . However, we can use 
X-ray colors to separate the points sources into g roups that are 
likely dominated by c ertain source types (e.g., IColbert et al.l 
120041: iPrestwich et ani2003l 12009,). We cons t ructed an X- 
ray color-color plot following IPrestwich et al.l (l2003h . The 
X-ray colors are defined as HC = (H-M)/T (hard color) 
and SC=(M-S)/T (soft color), where S, M, H, and T are 
the count rates in the 0.3-1.0 keV, 1.0-2.0 keV, 2.0-8.0 keV, 
and 0.3-8.0 keV bands, respectively. A Bayesian approach 
was adopted to calculate the colors and the associ a ted er - 
rors following the method developed by iPark et al] (l2006l) . 
This approach provides a rigorous statistical treatment of the 
Poisson nature of the detected photons as well as the non- 
Gaussian nature of the error propagation. The color-color 
plot is shown in Figure [8] We consider a source to be sig- 
nificant in the plot (a black data point) if its average error of 
the soft and hard colors is smaller than the 3(T-clipped mean 
(0.18) of the errors for all the sources. The other sources 
(cyan data points) have relatively large color errors and are 
not included in the following analysis. A significant fraction 
(75% of the bulge sources and 65% of the ring sources) of 
the X-ray sources a re located in the regi on that is likely dom- 
inated by LMXBs (|Prestwich et al.P 2003). The dominance of 
LMXBs is also suggested by the low SFR of NGC 1291. The 
SFR of NGC 1291 is ^ 2 0% of that of the Milky Way (e.g., 
IChomiuk & Povichll20Tlh . and the Milky Way only hosts a 
few HMXBs with luminosities above our detection complete- 
ness limit (« 10^^ erg s~^; se e Section 5.1) give n its popu- 
lation studies (e.g..lGri mm et a l. 2002; Voss & Aie llo 2010). 
Therefore, HMXBs should have a small contribution to the 
observed XRB population in NGC 1291. 

The bulge and ring sources appear to have different distri- 
butions in the X-ray color-color space. A two-dimensional 
Kolmogorov-Smirnov test''* indicates that the probability of 
the two distributions being drawn from the same parent pop- 

We used IDL routine KS2d for this test; see 



Fig. 8. — X-ray color-color plot for the point sources. 5, M, H, and T are 
the count rates in the 0.3-1.0 keV, 1.0-2.0 keV, 2.0-8.0 keV, and 0.3-8.0 keV 
bands, respectively. Black data points represent sources with relatively small 
color errors, while cyan data points represent sources with large eiTors and 
are less significant. The green and red tracks show the expected X-ray colors 
for absorbed power-law spectra with different power-law indices and col- 
umn dens ities. The blue e llipse indicates the area that is likely dominated by 
LMXBs iPrestwich et alj r2003). A significant fraction of the X-ray sources 
with small color errors are located in the region expected to be dominated by 
LMXBs. (A color version of this figure is available in the onhne journal.) 

ulation is 0.08. The 3(T-clipped mean soft and hard colors for 
the bulge population are 0.04 ±0.02 and -0.14 ±0.02, and 
the mean values are 0. 1 1 ± 0.02 and -0.15 ± 0.02 for the ring 
population. On average, the ring population appears to have 
a significantly harder soft X-ray color than the bulge popula- 
tion. The color-color plot is overlaid with a grid showing the 
expected locations of sources with different power-law spec- 
tra and absorption column densities, derived using the CXC's 
Portable, Interactive, Multi-Mission Simulator (PIMMS). A 
generally harder soft color suggests that the ring sources are 
more obscured (larger A^h) than the bulge sources, as would 
be expected in regions with ongoing/recent star formation. 

5. X-RAY LUMINOSITY FUNCTIONS 

We construct full-band XLFs for the bulge and ring regions 
using the merged Chandra observation. Although « 25^0% 
of the point sources are variable between the individual ob- 
servations, the XLFs derived from the merged observation 
should be stable against source variability base d on previous 
studie s of the X LFs for various systems (e.g., iGrimm et alj 
20051 I ZezaseTa l. 2007; Seh et al. 2011). 

5.1. Incompleteness Corrections and Background AGN 
Subtraction 

To derive the XLFs correctly, it is essential to make proper 
corrections for the detection incompleteness and background 
AGN contaminat i on. W e performed simulations following 
iKim & Fabbiand (1200 4*) to correct for the effects of incom- 
pleteness, including source detection limit and source confu- 
sion. Briefly, in each simulation we added a mock source at 
a random location of the event file of every observation using 
the MARX ray-tracing simulator,'^ created a merged event 
file and a merged full-band image, and then ran WAVDETECT 
to test if this additional source is detectable. The input X-ray 
luminosity of the source was randomly drawn from a power- 
law XLF with /3 = 1 in a cumulative form, A^(> Lx) = kL^ . 

|http://www.astro. Wash ington. edu/users/yoac him/idl/py_idl.html| 
http://space.mit.edu/CXC/MARX/index.html 
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We assumed a power-law spectrum for the source with F = 
1.7. The position of the source was randomly determined fol- 
lowing the r^l^ law ("de Vauco uleursl[T94 8'). We note that the 
adopted luminosity and position distributions here do not af- 
fect the incompleteness corrections significantly, as we only 
aimed to derive the positional-dependent detection fractions 
at a given luminosity. We performed 60,000 simulations in 
total, and the incompleteness corrections were calculated as 
a function of X-ray luminosity for the bulge and ring re- 
gions, respectively, utilizing all the simulated sources in each 
region. We also excluded the central 10"-radius area from 
the bulge region when computing XLFs as the incomplete- 
ness corrections are not reliable in this crowded region. The 
resulting 90% (50%) completeness limits are ss 1.5 x lO-*^ 
(w 7.1 X lO^^*") and « 2.2 x 10^^ (w 8.5 x lO^^*") erg s"' for 
the bulge and ring populations, respectively. The complete- 
ness limits in the bulge are lower because of the longer effec- 
tive exposure time in the bulge and the smaller background 
within the extraction regions of bulge sources. The XLFs 
were only calculated down to the 50% completeness lim- 
its, below which the Eddington bias (E ddington 1913) could 
have affected the measured X-ray luminosities. There are 64 
and 61 sources with luminosities above these completeness 
limits in the bulge and ring regions, respectively; w 98% of 
the sources have more than 10 full -band counts and thus the 
Eddington bias should have minimal effec t (e.g., Bauer et al] 
120041). Compared to the X LFs presented in llrwin e t al. (2002) 
and i Kilgard et al.l (l2005h . which were based on the first two 
Chandra observations, we are probing the XLFs about twice 
as deep. 

To correct for background AGN contamination, we did not 
attempt to classify ACiNs from the detected sources. Instead, 
we calculated the expected AGN luminosity function, which 
was then deducted from the observ ed XLFs. \ye co mputed the 
AGNs flux distribution using the lGilli et al.l (^2007) AGN PS 
model, and we normalized the AGN surface density to the ob- 
served value in the « 4 Ms Chandra Deep Field-South, which 
is ft^ 10000 deg" ^ at a flux Umit of w 3.2 x 10"'^ erg cm'^ s"' 
dXue et alj201 ll) . A 20% uncertainty was included for this es- 
timation that was attributed to cosmic variance (e.g. jLuo et alJ 
|2008). There are 4.2 ±0.8 and 21.5 ±4.3 AGNs expected 
in the bulge and ring regions, respectively. The ring XRB 
population is heavily contaminated (« 35%) by background 
AGNs. To derive the differential XLFs of the point sources 
in NGC 1291, we calculated the expected number of back- 
ground AGNs in a given luminosity bin and then subtracted it 
from the incompleteness-corrected number of X-ray sources 
in this bin. The cumulative XLFs were constructed by inte- 
grating the differential XLFs. 

The incompleteness-corrected and background subtracted 
XLFs for the bulge and ring populations are displayed in Fig- 
ure|9] The Icr errors were determined by combining the Pois- 
son errors of the number of sources in each luminosity bin and 
the 20% uncertainty of the AGN surface density. 

5.2. Parameterization of the X-Ray Luminosity Functions 

The shapes of XLFs for different XRB populations (mainly 
LMXBs and HM XBs) have been studied extensively in the 
literature (see e.g.. lFabbianol2006l for a review). For LMXBs, 
the XLF in the bright end (Lx > a few 10^^ erg s"') can 
usually be parameterized with broken power laws. The 
break luminosity is « (2-5) x lO-'^ erg s"', which may 
be related t o the Eddington luminosity of neutron star bi- 
naries (e.g.. ISarazin et alj|2000t iFinoguenov & JonesI 120021 : 



iKim & Fabbiand |2004|) . Below the break luminosity, the 
XLFs have differentia l power-law slopes of 1.8-2.2 (e.g., 
iKim & Fabbiand 12004 '): above the break, the XLFs cut off 
sharply, with differential slopes of > 2.8. Note that the 
XLFs for young (< 5 Gyr) LMXBs samples may not present 
this high-luminosity break (,Kim & Fabbiano 20 10). In the 
faint end (Lx < a few 10^^ erg s"^), there are limited data 
available for the XLFs of LMXBs. The XLFs of neai'by 
early-type galaxies and the bulges of nearby spirals (includ- 
ing the Milky Way) suggest that there is a second break at 
w (2-5) X 10^^ erg s~', below which the XLFs appear flat- 
ter with differential slopes of ~ 0-5-1 (e.g., Gilfanov 20041; 
IVoss & GilfanovllIOOTtlKrm et alJl2009l:IZhang et aL.201 1,) .'^ 
The origin of this low-luminosity break is not clear, which 
could be relat ed to magnetic braking during binary stellar 
evolution (e.g. JPostnov & Kuran ov 2005; Fragos et al. 2008^ 
or differe nt types of donor stars in high- and low-luminosity 
LMXBs dRevnivtsev et alj 1201 Ih . For HMXBs, the XLFs 
can be described with a simple power law with a differen- 
tial slope of w 1 .6 over a lar ge luminosity range (« 5 x 
1036_io4" erg s-^ lGrimm et"ar. 2003). 

We fit the XLFs of NGC 1291 to parameterize their shapes. 
A single power law of the form 



dN 



--A1 



^10" 



and a broken power law of the form 



dN 
dLx 



Al 
Al 



ifLx<Lb 
if Lx > Lb 



(2) 



(3) 



were used, where Al and Al are the numbers of sources at 
the reference luminosites. We used the SHERPA spectral fit- 
ting tool in the CIAO package to fit the XLFs in their differ- 
ential forms. It has been shown in previous studies that this 
is an efficient method to derive the characteristics of XLFs 
and it can yield consistent results with a maximurn likelihood 
appro ach (e.g., IZezas et al.ir2007l: iKim et al.1 120091: ISell et al.1 
1201 Ih . The XLFs were read into SHERPA as typical X-ray 
spectra with the number counts binned in small luminosity 
bins, (51og(Lx) = 0.05. Given the small numbers of sources 
in the bulge and ring regions, we used the Cash statistic in 
SHERPA for the model fits. The incompleteness corrections 
were supplied as an ancillary response file (ARE), and the 
AGN contribution was supplied as a table model in addi- 
tion to the intrinsic power-law model. The best-fit parameters 
with 1 a errors are reported in Table 21 and the best-fit bro- 
ken power-law models are shown in Figure |9l For the single- 
power model, both XLFs have a power slope of « 1.5. For 
the broken power-law model, both XLFs show a break lumi- 
nosity of w 5 X lO-'^ erg s"' . Above the break luminosity, the 
power slope is « 2; below that, the bulge XLF has a slope of 
w 1 while the ring XLF appears flatter, with a slope of « 0.4. 
Note that the errors for the broken power-law parameters are 
relatively large. 

The best-fit parameters were derived using the Cash statis- 
tic, which does not provide a goodness-of-fit parameter an- 
alytically. In order to compare the quality of the power- 
law and broken power-law fits, we adopted the likelihood 

It has been suggested that LMXBs in different environments (e.g., in the 
field or globular clusters) probably exhibit notable difference in their XLFs 
(e.g. . [Kim et alj.20Q^) . Such detailed discussion is beyond the scope of the 
cun'ent paper. 
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Fig. 9. — Differential [panel (a)] and cumulative [panel (b)] XLFs for the bulge (red filled dots or curve) and ring (blue open dots or curve) populations. The 
XLFs have been corrected for detection incompleteness and AGN contamination. In panel (a), the dashed curves represent the expected XLFs of background 
AGNs for the corresponding regions. The ring XRB population is heavily contaminated by background AGNs, resulting in large error bars for the XLF data 
points. The solid curves are the best broken power-law fits to the con'esponding XLFs. In panel (b), the red dotted curves represent the Icr errors for the ring 
XLF, and the grey region represents the Icr errors for the bulge XLF. The cumulative XLFs do not appear to monotonically decrease with luminosity at a few 
places due to the uncertainties introduced by the background subtraction. The XLFs were computed down to the 50% completeness limits. (A color version of 
this figure is available in the online journal.) 



ratio test used in IZezas et all (12007b . For each differential 
XLF (ring/bulge), we simulated 10,000 mock XLFs from the 
best-fit power-law model, using the same luminosity bins and 
adding Poisson fluctuations to the source counts. The de- 
tection incompleteness and background AGN contamination 
were also applied to the mock XLFs. We then fit the simulated 
XLFs with both the power-law and broken power-law models. 
The ratio of the Cash statistic values of the best power-law 
and broken power-law fits for the simulated XLFs were cal- 
culated and compared to the ratio for the observed XLF; the 
rate of obtaining a higher ratio in the simulations represents 
the probability of getting an improved broken power-law fit by 
chance. For the bulge XLF, we found that the broken power- 
law model provides a statistically improved fit at the 89% con- 
fidence level; for the ring XLF the broken power-law model 
provides an improved fit at the 94% confidence level. There- 
fore, the broken power-law model appears to provide better 
fits to the data for both the bulge and ring populations. The 
broken power-law shapes of the XLFs suggest that both the 
bulge and ring XRB populations are dominated by LMXBs. 
We are probing the low-luminosity break at w 5 x lO"*^ erg s"' , 
with a slope of w 2 at the brighter end and a flatter shape at the 
fainter end; there are too few luminous XRBs in NGC 1291 to 
constrain the high-luminosity break at w (2-5) x lO-'** erg s"' . 
This result is consistent with that inferred from the X-ray col- 
ors (see Section 14.41) . Given the arguments above, we pre- 
fer the broken power-law model for both the bulge and ring 
XLFs, and the following discussions are based on the broken 
power-law parameters. 

6. THEORETICAL MODELING OF THE X-RAY BINARY 
POPULATIONS 

PS modeling has been proven a valuable tool in order to 
study the physical properties of the XRB populations ob- 
served with Chandra in nearby galaxies. Belczynski et alJ 
(l2004l) constructed the first synthetic XRB populations for 
direct comparison with the Chandra observed XLF of NGC 
1569, a_star forming dwarf irregular galaxy. iLinden et alJ 
(|20()9, 2010) developed models for the HMXB and Be XRBs 
of the SMC, studying the XLF and the spatial distribution of 
this population and investigating the effect of electron-capture 



supernovae of massive ONeMg stellar cores. iFragos et al.l 
( l2008i 12009 ) performed extensive PS simulations for the 
modeling of the two old elliptical galaxies NGC 3379 and 
NGC 4278. They studied the relative contribution to the ob- 
served XLF from sub-populations of LMXBs with different 
donor and accretor types, and the effects of the transient be- 
havior of LMXBs. All these studies have been pivotal in in- 
terpreting Chandra observations of nearby galaxies. 

In this paper, we utilized a new grid of PS models to in- 
terpret the XLF characteristics of NGC 1291. In these mod- 
els, we focused on LMXBs and HMXBs formed in the galac- 
tic field as products of the evolution of isolated primordial 
binaries (i.e. we did not take into account any dynamical 
effects). We performed our simulations with StarTrack 
(iBelczynski et al.ll2002L l2008h . an advanced PS code that has 
been tested and calibrated using detailed mass transfer star 
calculations and observations of binary populations, and in- 
corporates all the important physical processes of binary evo- 
lution. 

In the development of our models we incorporate our cur- 
rent knowledge about the characteristics of the stellar pop- 
ulation in NGC 1291. This galaxy has a stellar mass of 

1.5 X 10" MfT, and a low current SFR of 0.4 yr~' 
(iKennicuttetaLl 120031) . We adopted the iNoU et all (l2009l) 
model for the star-formation history of NGC 1291 including 
two exponential starburst events (see Section 1). The metallic- 
ity of both populations is assumed to be solar. Due to the sim- 
ilar shapes of the bulge and ring cumulative XLFs and their 
relatively large uncertainties, it is not plausible to compare 
models to the individual XL Fs. We also note that the best fit 
results by iNoll et af] (l2009l) about the star-formation history 
of NGC 1291 have no spatial information. Hence, our com- 
parison between PS models and the observed XRB population 
is limited to the combined (bulge plus ring) XLF of the galaxy. 

We created a grid of 48 PS mo dels varying pararn eters that 
are known from previous studies dFragos et al.ll2008l) to affect 
the shape and normalization of the resulting synthetic XLF. 
Namely, we varied the common envelope (CE) efficiency, the 
initial mass function, stellar wind strength, and finally we al- 
lowed for various outcomes of the CE phase, in particular the 
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possible CE inspirals with Hertzsprang gap donors that ter- 
minat e binary evolution bar ring the subsequent XRB forma- 
tion ( Belczvnski et al.|[2007h . The latter parameter, which is 
known to affect the evolution of massive bin aries which are 
the pr ogenitors of double compact objects (Belczvns ki et alJ 
l200 7h. is for the first time explored in the context of XRB 
populations. The various parameter values used are shown in 
table|5] We note that in our calculations we combined acE and 
A into one CE parameter, where A is a measure of the central 
concentration of the donor and the envelope binding energy. 
In the rest of the text, whenever we mention the CE efficiency 
acE, we refer in practice to the p roduct acE x A, effectiv ely 
treating A as a free parameter (see lBelczynski et"ani2008l for 
details). 

The synthe tic XLFs wer e constructed based on the method 
described in iFrag os et al.l (1200 8 ), wh ile transient systems 
were modeled according to lFragos et al.i (2009). Furthermore, 
we calculated the errors in the synthetic XLF due to Poisson 
low number statistics (high luminosity end) and uncertain- 
ties in the bolometric correction factors (low luminosity end) 
(T. Fragos et al. 2012, in prep.; R Tzanavaris et al. 2012, in 
prep.). In order to assess which are the preferred values of the 
different model parameters, we calculated the likelihood of 
the observed XLF given a model. Our statistical analysis takes 
into account not only the shape of the XLF (as in Fragos et al.l 
HbOS) but the absolute normalization too. Our maximum like- 
lihood model has a low CE efficien cy (acE ^ 01) and the 
"standard" stellar wind prescription (iBelczvnski et al.l [20101 
and references therein). We also found that there is an appar- 
ent degeneracy between the slope of IMF and whether or not 
we restrict the outcome of a CE event to a merger when the 
donor star is in the Hertzsprung gap. A flatter IMF results in 
a larger black hole XRB (BH-XRB) population, while on the 
other hand, assuming that a CE with a Hertzsprung gap donor 
leads always to a merger, which prevents many of the binaries 
from evolving into a BH-XRB. Hence, a BH-XRB population 
consistent with observations could either result from a steeper 
IMF (slope of -2.7) and CE prescription where we allow for 
all possible outcomes, even with a Hertzsprung gap donor, 
or it could result from a flatter IMF (slope of -2.35) and a 
CE model where a CE event with a Hertzsprung gap donor 
would always lead to a merger. In fact, the former case corre- 
sponds to our maximum likelihood model, while the latter to 
the model with the second highest likelihood. 

Figure [TO] shows the synthetic XLF from our maximum 
likelihood model in comparison with the observed XLF, as 
well as the different types of XRBs contributing to the XLF 
at different ranges of X-ray luminosity. There are two ar- 
tificial "jumps" in the model XLF at w 10^^ erg s"' and 
« 8 X 10^'^ erg s"', instead of the observed smoother power- 
law shape. In our models, we identify transient and persistent 
sources by comparing the calculated mass-transfer rate of the 
binary to a critical mass-transfer rate, below which the ther- 
mal instability develops. Furthermore, we strictly limit the 
accretion rate to the Eddington limit. These limiting mass- 
transfer rates in our modeling are responsible for the apparent 
jumps in the synthetic XLF. However, in nature there is no 
sharp transition between thermally stable and unstable disks, 
nor a precise limit in the highest accretion rate possible. Ac- 
cretion onto a compact object is a non-linear and much more 
complex process, which in reality can results in a smoother 
luminosity distribution with no sharp transitions. Finally, by 
assuming a solar metallicity for the s tellar popu l ation o f NGC 
1291 in order to be consistent with iNoll et al.l (l2009l) and in 



the absence of any stellar metallicity measurements, we im- 
pose in practice a maximum BH mass of « 15 Mq. If instead 
part of the stellar population had a lower metallicity (e.g. 30% 
solar) th en the maximum BH m ass would increase to 30 M© 
or more (iBelczynski et al.ir2010 ). which could smooth out the 
jump at the very luminous end of the XLF. 

We found that the BH-XRBs that populate the high-end of 
the XLF (above « 2 x 10^** erg s"', Eddington limit for a 
hydrogen accreting neutron star) are transient systems with 
evolved companion stars, while at lower luminosities the 
XRB population consists of a mixture of transient and per- 
sistent XRBs with a neutron star accretor and a giant or a 
main sequence (MS) donor star. These findings and the gen- 
eral shap e of the synthetic X LF are in agreement with the PS 
study by Fragos et al. (2008) of the two old elliptical galax- 
ies NGC 3379 and NGC 4278. This leads us to the con- 
clusion that an old LMXB population is dominating the ob- 
served XLF. However, there is an evident overabundance of 
super-Eddington (Lx ^ 2 x 10^^ erg s"', the Eddington limit 
for a hydrog en accreting 1 .4 M© neutron star) XRBs, com- 
pared to the iFragos et al.l (l2008l) models, which can be at- 
tributed to the enrichment of the younger XRB population 
with black hole accretors. This younger population is likely 
originated from the ring star-forming region. Our finding that 
the younger stellar population of the ring region enhances 
the high X-ray luminosity, super-Eddington BH-XRB pop- 
ulation, is consis t ent w it h the recent obser vational study by 
iKim & Fabbiand (|2010|) . iKim & Fabbianol (|2010) compared 
the combined XLF of LMXBs detected in Chandra observa- 
tions of post-merger elliptical galaxies with that of typical old 
elliptical galaxies and found that elliptical galaxies with signs 
of recent star formation host a larger fraction of luminous 
LMXBs (Lx > 5 X 10^^ erg s"') than old elliptical galaxies. 

7. DISCUSSION 

7.1. The Bulge LMXB Population 

Galactic bulges host relatively homogeneous samples of 
old stellar populations, and thus they are ideal targets for 
the research of LMXB XLFs. PS modeling of the XLFs has 
been used to constrain the formation, evolution, and physical 
prope rties of the LMXB population (e.g., Fragos et al. 20081 
|2009|) . Recently deep Chandra exposures of nearby ellipti- 
cal galaxies and spiral bulges have pushed the XLFs down 
to limiting luminosities of around a few 1 0^^ erg s~^; some 
best studied examples include Centaurus A (Voss & Gilfanovl 
2006), M31 ( Voss & Gilfanov 2007), M81 (Sell et al. 2011), 
and samp l es of galactic bulges in iGilfanovl {2004) and 
iKim et al] (120091) . The derived XLFs are remarkably simi- 
lar, with a low-luminosity break at w 5 x lO^*^ erg s"' and a 
high-luminosity break or cut-off at « 5 x 10^^ erg s"'. The 
bulge of NGC 1291 is dominated by LMXBs based on the 
XLF shape and the X-ray colors, though relatively weak star 
formation could still exist given the observed filaments of Ha 
gas (See Section 1). The bulge XLF of NGC 1291 reaches a 
similar threshold luminosity of 7 x 10^'' erg s"', and its shape 
is statistically consistent with previous findings. We obtain a 
low-luminosity break at « 5 x lO-'^ erg s"'; a high-luminosity 
break is not constrained but is highly probable as there is only 
one bulge source with > 5 x 10^^ erg s"' (one of the ULX 
candidates). The power-law slope above the low-luminosity 
break is w 2, and is ss 1 below the break, typical for LMXB 
XLFs. Despite of the complex morphology of the galaxy, the 
bulge of NGC 1291 stifl hosts a population of LMXBs that 
is very similar to those in typical elliptical galaxies and spiral 




Fig. 10. — Left panel: Synthetic XLF (black solid curve) for the maximum likelihood model. The gray shaded area corresponds to the l—cr error of the XLF 
due to Poisson statistics and uncertainties in bolometric correction. For comparison the observed combined (corresponding to both the bulge and ring regions) 
XLF (red solid curve) and the associated l—cr error (red dotted curves) are drawn. Center panel: Contribution of different sub-populations by separating the 
XRBs into groups of systems with different donor stellar types. Right panel: Contribution of different sub-populations by separating the XRBs into groups with 
different accretor types. In the center and left panel, solid colored lines denote persistent XRB sub-populations, while dashed colored lines denote transient ones. 
(A color version of this figure is available in the online journal.) 

in Figure [TT^. The numbers of sources have been cor- 
rected for the detection incompleteness and background 
AGN contamination. The fraction of relatively lumi- 

erg s" 



bulges (in terms of their XLFs) down to the limiting luminos- 
ity, suggesting that all LMXB populations are compatible. 

7.2. Comparison of the Bulge and Ring Populations 

X-ray observations of the XRBs can provide valuable con- 
straints on the underlying stellar populations, independent of 
what were obtained from the other wavelengths. Given the 
X-ray observations of NGC 1291, we consider that the ring of 
this galaxy hosts a younger stellar population than the bulge, 
based on the arguments given below. 



1. 



There appears to be an overdensity of X-ray sources 
in the ring of the galaxy. The radial profile of the X- 
ray source surface density follows the NUV surface 
brightness distribution in the ring instead of the /T-band 
surface brightness distribution, suggesting that the ring 
XRBs are likely associated with a younger stellar pop- 
ulation (see Section l4n) . 



2. The ring XRB population has in general a harder soft 
X-ray color than the bulge population, likely due to in- 
creased obscuration in the ring associated with ongo- 
ing/recent star-forming activities (see Section l4!4b . 

3. The bulge and ring combined XLF can be well mod- 
eled with a synthetic XLF that is in general agreement 
with the PS modeling of typical old elliptical galaxies. 
However, there appears to be an enhancement of super- 
Eddington XRBs, which could be associated a possible 
younger stellar population in the ring region (see Sec- 
tion|6l). 

4. The XLF is an efficient probe of the relative age of stel- 
lar populations. It has been reported that young ellipti- 
cal galaxies host a larger fraction of luminous LMXBs 
and the XLFs app ear flatter than th ose of typical old 
elliptical galaxies dKim & Fabbiand 12010) . For NGC 
1291, the slope of the broken power-law fit at the faint 
end (TablelDl for the ring appears flatter than that for the 
bulge, albeit with a large uncertainty, suggesting that 
there are more luminous X-ray sources in the ring. A 
comparison of the numbers of sources in given lumi- 
nosity bins for the bulge and ring populations is shown 



nous sources (Lx > 5.5 x lO"*^ erg s ')) is 30^% for 
the bulge, and it is 54;!;[|% for the ring; the errors are 
Poissonian. 

The LMXBs of the ring likely belong to two popula- 
tions, a relatively old population that ages similarly to 
the bulge population and a relatively young population 
that formed with the enhanced star formation. To iso- 
late this young population, we simply assume that the 
source number counts of the old population scale with 
the total A'-band luminosity, and then the old popula- 
tion in the ring can be derived from the bulge source 
number counts and the ratio between the /T-band lumi- 
nosities of the bulge and ring (a factor of 3.6). We sub- 
tract this old population from the ring number counts 
and the residual is displayed in Figure [TTb. The frac- 
tion of relatively luminous sources is 12^^% for this 
remaining population in the ring. The difference in the 
fraction of luminous sources between the bulge and the 
ring appears more significant, and it is likely that the 
remaining XRB population is indeed younger. 

7.3. XRB Population of Galactic Rings 

As the two best studied galactic rings in the X-ray, the 
rings of NGC 1291 and the Cartwheel galaxy differ signif- 
icantly. The XRB population in the ring of the Cartwheel 
galaxy is dominated by ULXs; the XLF has a power-law 
slope of sa 1.6 albei t at high luminosities (Lx 10"''^ erg s"' ; 
IWolter & Trinchieril|2004j) . consistent with being dominated 
by luminous young HMXBs. For NGC 1291, there are only 
two ULXs in the ring, and the ring XLF shape and X-ray 
source colors suggest a LMXB dominated system. The ring 
of the Cartwheel galaxy clearly has much stronger ongo- 
ing star formation than the ring of NGC 1291 (total SFR of 
0.4 Mq yr"' for NGC 1291 compared to 67.5 M© yr"' for the 
Cartwheel galaxy). 

Such a discrepancy is probably caused by different forma- 
tion and evolution scenarios of the galactic rings. The ring of 
the Cartwheel galaxy is likely produced by a recent merger 
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Fig. 1 1 . — Numbers of X-ray sources in given luminosity bins for the bulge, ring [panel (a)], and bulge-subtracted ring [panel (b); see text] populations. The 
numbers of sources have been corrected for the detection incompleteness and background AGN contamination. The luminosity bins were chosen to be the same 
for all the populations, and the limiting luminosity is Ri 8.5 X 10^' erg s"' . The luminosities for the bulge data points were shifted by a small amount for display 
purpose. For several cases with only error bars and no data points shown in the plot, the numbers of sources are smaller than 0.1. The vertical dashed lines mark 
the approximate breaks of the XLFs, and are used to separate relatively luminous and relatively faint sources. (A color version of this figure is available in the 
online journal.) 



event, and the expanding density wave forms the starburst 
ring. CoUisional rings are generally short-lived with a lifetime 
of the order of a few hund red million years (see the review of 
iMoiseev & Bizyaevll2009l and references therein). The ring of 
NGC 1291 is probably created by secular gas accumulation 
at the resonance orbit due to the gravity torques of the bars. 
Enhanced star formation appeared in the ring as gas was be- 
ing gathered. However, this could a be a less energetic and 
longer-term event than those appeared in collisional ring sys- 
tems, resulting in th e relatively mild s tar formation observed 
today. Adopting the Noll et all (12009*) star-formation history 
model (Section 1) and assuming the ring hosts the young stel- 
lar population with an e-folding timescale of 50 Myr, the SFR 
of NGC 1291 would be w 20 Mq yr"' 200 Myr ago, compa- 
rable to the current SFR of the Cartwheel galaxy. 

8. SUMMARY 

We presented Chandra observations of the ring galaxy 
NGC 1291 and studied the XRB populations in the ring and 
bulge regions of the galaxy. The main results are summarized 
below. 

1. We analyzed the four Chandra AC1S-S3 observations 
of NGC 1291, with a total exposure of w 180 ks. We 
have detected 169 point sources in the 0.3-8.0 keV 
band, 75 of which are in the bulge and 71 in the ring. 
The full-band luminosities of the sources range from 
2.7 X 10^'' to 2.3 X 10^''ergs-i. 

2. We studied basic properties of the X-ray sources. There 
appears to be an overdensity of X-ray sources in the 
ring of the galaxy. The radial profile of the X-ray source 
surface density resembles the NUV surface brightness 
distribution in the ring, suggesting that a significant 
fraction of the ring XRBs are associated with a rel- 
atively young stellar population. About 40% of the 
sources in the bulge and 25% of the sources in the ring 
show > 3ct long-term variabiUty in the full-band count 
rate. 

3. We studied the spectra of the nuclear source, which ap- 
pears to be a low-luminosity AGN with moderate ob- 
scuration, and there is also hot gas present in the nuclear 



region that produces soft excesses in the spectra. The 
spectra of individual observations were well fit with an 
absorbed power-law plus thermal plasma model. 

4. In the X-ray color-color plot, a significant fraction of 
the bulge (« 75%) and ring (« 65%) sources with rel- 
atively small color errors are located in the region ex- 
pected for LMXBs. On average, the ring XRB popu- 
lation has a harder soft X-ray color, suggesting that the 
ring XRBs are more obscured, probably associated with 
ongoing/recent star-forming activities in the ring. 

5. We constructed full-band XLFs for the bulge and ring 
regions, taking into account the detection incomplete- 
ness and background AGN contamination. The re- 
sulting 90% completeness limits are w 1.5 x lO-'^ and 
« 2.2 X 10^' erg s"' for the bulge and ring populations, 
respectively. The ring XRB population is heavily con- 
taminated (sa 35%) by background AGNs. Both XLFs 
can be fit with a broken power-law model. The shapes 
are consistent with those expected for populations dom- 
inated by LMXBs. The ring XLF appears flatter than 
the bulge XLF, indicating that the ring hosts a larger 
fraction of luminous sources and is likely associated 
with a younger stellar population. 

6. We successfully modeled the bulge and ring combined 
XLF with a PS model, utilizing our current knowledge 
of the galaxy properties. The model is consistent with 
those for typical old elliptical galaxies, suggesting that 
an old LMXB population is dominating the observed 
XLF. The relative enrichment of luminous sources of 
NGC 1291 can be attributed a younger XRB population 
originated from the ring star-forming region. 

7. We consider that the ring of NGC 1291 hosts a younger 
stellar population than the bulge based on the X-ray 
studies presented in this paper. 
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TABLE 1 

Chandra Observations OF NGC 1291 



Obs. ID 


Start Date 


Exp (ks) 


Cleaned Exp (ks) 


PI 


(1) 


(2) 


(3) 


(4) 


(5) 


795 


2000 Jun 27 


39.7 


37.6 


J. A. Irwin 


2059 


2000 Dec 7 


37.0 


22.8 


A. H. Prestwich 


11272 


2010 May 4 


70.0 


69.1 


G. Fabbiano 


11273 


2010 May 8 


50.0 


49.1 


G. Fabbiano 



Note. — Column 1: Chandra observation identification number Column 
3: Observation start data. Column 3: Nominal exposure time. Column 4: 
Exposure time after removing background flares. Column 5: Name of the 
Principal Investigator. 



TABLE 2 
Chandra Point-Source Catalog 



X-ray Coordinates 



Counts 



XID 


0:2000 


52000 


Off-Axis 




FB 


FB Low Err 


FB Up Err 


SB 


SB Low Err 


SB Up Err 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(V) 


(8) 


(9) 


(10) 


(11) 


1 ... 


03 16 47.17 


-41 08 15.2 


4.6 


6.2 


29.9 


5.9 


7.2 


19.2 


4.8 


6.2 


2... 


03 16 47.62 


-41 05 34.2 


4.2 


5.9 


144.8 


13.0 


14.2 


106.9 


11.3 


12.6 


3 ... 


03 16 48.09 


-Al 04 10.1 


4.6 


6.2 


6.6 


2.8 


4.2 


10.8 


-1.0 


-1.0 



Note. — Units of right ascension are hours, minutes, and seconds, and units of dechnation are degrees, arcminutes, and arcseconds. (This table is 
available in its entirety in a machine-readable form in the onhne journal. A portion is shown here for guidance regarding its form and content. The full 
table contains 34 columns of information on the 169 X-ray sources.) 



TABLE 3 

Best-Fit Parameters for the Nuclear Source 



Observation x^/dof 



(1) 



(2) 



r A'h,2 ^-2-8 

(X 10^^ cm-2) (X 10^' erg s"') 
(3) (4) (5) 



795 
2059 
11272 
11273 



0.84 
0.93 
1.21 
1.11 



1.6: 
2.6 

2.1! 
1.8 



2 0+*" 

i.8!Si 



'-0.5 



1. 



'-0.6 



2 57+0.08 
1 00+8:^^ 

I./S_g51 

1 90+0.04 
-0.23 



1.27 



NOTIi. — The spectra were tit with an absorbed power law plus thermal plasma 
model (WABSl*APliC+WABS2*P0W). N\ip_ and F are the absorption column den- 
sity and photon index for the power-law component, respectively. The temperature 
(0.15 keV) and column density (A'h.i =7.6x 10^^ cm""^) of the thermal component 
were fixed at the values derived from the stacked spectrum. The observed 2.0-8.0 
keV luminosities are dominated by the power-law component. 



TABLE 4 

Best-Fit Parameters for the X-Ray Luminosity Functions 



Region Single Power-law 



Broken Power-law 



Al 



a2 



Bulge 
Ring 



143 
159 



1 6-*<'i 
-Ql 
1 4+0.1 
^•^-0.1 



36!^ 
182 



138 
153 



1 l-**-^ 
0.4^1 



2 O+O-' 
^•^-0.3 



4.9«:^ 
6.4«.; 



26!j 
fi+T 



Note. — 5p and 5bp are the Cash fit statistic values for the power-law and broken 
power-law fit, respectively. See Equations (2) and (3) for the definition of the other 
parameters. 
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TABLES 
Model Parameters forNGC1291 



Pai'ameter 



Notation 



Value 



Star Formation History 

Metallicity Z 

Total Stellar Mass M, 

Binary Fraction F(,i„ 
IMF power-law slope 

CE Efficiency ocE 

Stellar wind strength fjwind 
CE during HG 

Magnetic Braking 



Noll et al. (2009) 
0.02 
1.5 X 10" M0 

50% 
-2.7 or -2.35 
0.1,0.2, 0.3, or 0.5 
0.25, 1.0, or 2.0 
Yes or No 
Ivanova & Taam (2003) 



